The need to smoothly cover a computational domain of interest genericallyrequires the adoption of several grids. To solve the problem of interest underthis grid-structure one must ensure the suitable transfer of information amongthe different grids involved. In this work we discuss a technique that allowsone to construct finite difference schemes of arbitrary high order which areguaranteed to satisfy linear numerical and strict stability. The techniquerelies on the use of difference operators satisfying summation by parts and{\it penalty techniques} to transfer information between the grids. This allowsthe derivation of semidiscrete energy estimates for problems admitting suchestimates at the continuum. We analyze several aspects of this technique whenused in conjuction with high order schemes and illustrate its use in one, twoand three dimensional numerical relativity model problems with non-trivialtopologies, including truly spherical black hole excision.
展开▼
